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We argue that recent data on fluctuations observed in heavy ion colli- 
sions show non-monotonic behaviour as function of number of participants 
(or "wounded nucleons") N\y. When interpreted in thermodynamical ap- 
proach this result can be associated with a possible structure occuring in 
the corresponding equation of state (EoS). This in turn could be further 
interpreted as due to the occurence of some characteristic points (like soft- 
est point or critical point) of EoS discussed in the literature and therefore 
be regarded as a possible signal of the QGP formation in such collisions. 
We show, however, that the actual situation is still far from being clear 
and calls for more investigations of fluctuation phenomena in multiparticle 
production processes to be performed. 

PACS numbers: 24.60.Ky; 25.75.Nq; 05.20.-y 



1. Introduction 

The primary goal of all high energy heavy ion experiments is to find the 
phase transition between hadronic matter and quark-gluon plasma (QGP) 
as predicted by the lattice QCD 1 . It is natural then to look at experimental 
results from such perspective and to scrutinize them for the possible signals 
of QGP formation. In this work we would like to concentrate on one charac- 
teristic feature, which could signal such phase transition, namely it should 

For relevant references see, for example, proceedings on recent QM conference 
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be accompanied by enhanced fluctuations in some variables [21 IS1 lU El • For 
it turns out that, when measuring in final states only hadronic probes pro- 
duced with small transverse momenta, the only observables surviving the 
possible phase transition (and, in principle, actually depending on its form) 
is the observed fluctuation pattern in properly selected variables. Following 
[2] we shall be interested in the quantity being the ratio of fluctuations in 
entropy to fluctuations in energy. The reason is that, in thermodynamical 
approach to be pursued here, both quantities are well defined for any form of 
matter (confined, mixed and deconfined) both in early stage of collisions and 
in the final state of the hadronizing system under consideration. When pro- 
duced matter can be treated as an isolated system, energy is conserved and 
entropy is also expected to be conserved during the expansion and freezeout 
stages of the hadronization. Therefore simultaneous measurements of both 
quantities are likely to provide us with information on the equation of state 
(EoS) of the hadronizing system. This ratio can be written (see Appendix 
A for details) as 

„ iff /, a V 2 , , d ^V TdV 

R=- — ^7= 1 + r > where — - = — — . (1 

' 6e\ V l + ach ' Y dlnT V dT y ' 



Here a = dp/de = p/e characterizes EoS, p denotes the pressure and 
e = E/V energy density with V being the volume of our system and T 
its temperature (statistical equilibrium is assumed). The entropy of the 
system, S, should be connected with the multiplicity of the produced sec- 
ondaries, S ~ N and, similarly, energy released in the production process, 
E, should be connected with the measured sum of transverse momenta, 
E ~ J2iLiPT- When variations of the reaction volume (identified with the 
volume of system) can be neglected, i.e., for = 0, one gets formula pro- 
posed in |2] as the right quantity the energy dependence of which should be 
investigated experimentally 2 : 



mf \ e) (1 + a) 2 y ' 

E' 2 

Notice that R is directly sensitive to the equation of state (EoS) of hadroniz- 
ing matter because of occurrence of parameter a. In what follows we shall 
assume that 

5S 5N 5E 5J2pt 

T=aT and -W=(E^)- (3) 



For recent reviews on this matter see [Sj- 
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We identify throughout this work 5N and SJ2pt with fluctuations caused 
by all possible sources (as only such are accessible in experiment), i.e., also 
all our results concerning different estimations oip/s by means of eq. (J2J) 
include influences of all possible sources of fluctuations. In 2 R has been in- 
vestigated by using some specific statistical model of nuclear collisions. The 
predicted very characteristic shape of energy dependence of this quantity 
(non-monotonic, with a single "shark fin" -like maximum) was then pro- 
posed as a possible signal of QGP phase transition which should therefore 
be subjected to future experimental tests jS]. 

In this work we would like to bring attention to the fact that apparently 
similar kind of fluctuations (but in function of Nyy, which can be connected 
with the reaction volume, rather than energy, as is the case in |2J) have 
been already observed and reported in [H] and In what follows, using 
some minimal input, we shall in the next Section rewrite results of [SI E] 
in terms of variable R defined by © and (JHJ), clearly demonstrating that it 
has similar non-monotonic character as that obtained in |2j. It is therefore 
tempting to argue that these data show us either what has been called in 
the literature the softest point of the corresponding EoS ^Oj or what is 
named in other investigations the critical point of EoS However, at 

this stage one not only cannot distinguish between both possibilities but 
even one cannot take them very seriously, for when confronting these data 
with similar data obtained recently in ^3], which use different measure 
of fluctuations, one discoveres that no such feature are seen there 3 . Still, 
there are some other possibilities offered recently, as, for example, spinodial 
decomposition occuring at the hadronization stage ^2], which deserve also 
attention. 

In Section III we shall also analyse (with the same aim as above, i.e., 
looking for a possible signal of QGP) recent interferometric data, which pro- 
vide us information on the interaction volume JKj 4 . It turns out that they 
lead to similar results for EoS as data from [SJ |§] (albeit this time at differ- 
ent energy). The last Section contains our summary whereas Appendices 
contain derivation of our basic formulas. 



3 Ref. |1M] provides data for parameter $ (defined in |14p at energies 40, 80 and 158 
AGeV for the following centrality bins: - 5%, - 6.5%, 5 - 10%, 10 - 15% and 
15 - 20%. 

4 The volume V measured here is the freeze-out volume but in our discussion we identify 
it with the production volume. 
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2. Nonmonotonical dependences observed in data on fluctuations 

In [HJ the following measure of fluctuations has been presented as a 
function of the number of participants Nyy (or "wounded" nucleons) 5 : 



,-, ^data ^random p. y (Pt) (Pt) 

T = o ' TfT~\ ' 

* ^random \PT) 

It increases with N\y, reaches maximum at Nw ~ 200, and decreases for 
higher values of Nyy. Whereas in [7] such behaviour has been attributed 
to the peculiar feature of the hadronization model used, namely to the 
percolation of hadronizing strings produced in collision process 6 , here we 
shall connect it directly, in a similar way as in [2], to the behaviour of EoS 
of the matter produced at the early stage of the collision. This can be done 
by rewritting Ft in terms of R defined in eq. © . To do so let us first notice 
that, because 

Var&PT) _ Var(p T ) 1 Var(N) 
(Ept) 2 (Pt) 2 (N) + (AO 2 ' 1 j 

one can write R as 

Var(N) 



Var(J2PT) 1 + 



where 



1 Var(N) , Var(p T ) 

r = ^^NT and " = 1m^- (7) 

On the other hand, as shown in Appendix B, Ft can be expressed in terms 
of r and parameter p, 



F T = V 1 + 2r - 2py/r 2 + r - 1, (8) 

where p € [—1,1], the correlation coefficient between number of particles Af 
and J^Pt = Ya=iPt, is our minimal input mentioned above 7 . Using © it 

5 For proper definition of different types of averages and mean values used here see 

EH 

6 Without introducing notion of EoS or phase transitions, unless phenomenon of perco- 
lation itself is regarded as a kind of phase transition. But even then it would be phase 
transition between bigger and smaller number of strings only, not between hadronic 
matter and quark-gluon plasma. 

7 It appears because, whereas R depends on fluctuations in N and on fluctuations in 
'YliLiPT, so far in experiment one measures quantities like $ or Ft in which both 
type of fluctuations occur together. Parameter p describes therefore their mutual 
(unknown a priori) correlations. 



rww-appb printed on February 1, 2008 



■5 



is now straightforward to rewrite eq. (JHJ) in the following form: 



(1 + F T ) 2 - 1 



R 



1 



2 1-R 
from which our main result follows: 

p 2 + 2y(l +y)±p- vV + 4y(l + y) 



R 



2(1 +y) 2 



where 



1 r 



V= 2 



(1 + F T ) 2 - 1 



(9) 



(10) 



11 



In addition to the a priori unknown correlation coefficient p there is also 
some freedom in the choice of sign in (|1U|) . To fix both of them let us notice 
the following: 

(a) It is known that purely statistical or broader fluctuations encountered 
in all multiparticle production processes, i.e., the fact that Var(N) > 
(N), lead to the condition that 8 



Var(N) __ uj 



> i. 



(12) 



(6) Furthermore, analysis of recent CERES data |13| shows clearly, see 
Fig. |1] below, that parameter oj is very slowly varying function of 
energy and multiplicity (esentially u — 0.43). This practical constancy 
of fluctuations in pt means that fluctuations of energy E, which are 
given by fluctuations of J2pt, is not so much given by fluctuations in 
Pt but by fluctuations in multiplicity N, cf., eq.©. 

These two observations mean therefore that for the multiplicity distributions 
of the poissonian type and broader, R is limited from below: 

R > — !— > 0.699 (13) 

(for uj = 0.43). Using this limit together with eq. (|10j) one obtaines that 
also p 2 is limited from below, namely 

? > (14) 

UJ + 1 



See, for example, |18|. Multiplicity distributions broader than poissonian distribu- 
tions are also observed in hadronic collisions |19|. 
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Because, without any additional correlations between x% and N, the Yli=i x % 
increases with N, the positive correlation coefficient, p > 0, seems to be the 
only natural choice. In addition, results of PHENIX [5] show that Ft is very 
small (actually not exceeding Ft ~ 0.04). Neglecting therefore term with y 
in (|14jl one can estimate that p > 0.83 (in our calculations we have put it 
slightly above this limit using value p = 0.85). It means that R corresponds 
to a stronger than poissonian fluctuations. To be consistent with limitations 
on R imposed by f|13|) we shall choose to the solution with positive sign in 
(fTU|) (the other one would lead to unacceptable small values of R) . 



0.75 
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Fig. 1. Left panel: results of transforming PHENIX data from jHj (i.e., Ft versus 
Nyy) to R as function of the number of participant (or number of "wounded") 
nucleons, Nw by using ea. H10|) . As explained in the text the positive sign in 1|10|) 
has been chosen and the correlation coefficient set to p = 0.85 (see text for details). 
Right panel: R obtained from 4> as measured by NA49 |§] (for the same value of 
p)- 

We are now prepared to translate, by using eq. (fTU|) . the experimental 
knowledge of Ft and $ as a function of number of participants provided 
in [SI into similar dependence of R. The results are presented in Fig. 
n The only unknown feature there is the choice of parameter p. In fact 
for p =const different values result in essentially the same shape of R, only 
shifted accordingly paralelly up or dow. The case of p = p(N\y) would 
lead to some changes, however. Unfortunately at the moment there are 
no data available to estimate the functional form of p{Ny/). We expect, 
however, that our choice of p used here corresponds to a lower limit for a 
possible effect observed in Fig. ^ To substantiate this we present in Fig. [2 
parameter R as function of N\y calculated directly from definition i.e., 
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using directly measured information on Var(N)/ (N) and Var(pT) / '(pt) 2 
obtained by [S] 9 . 



0.85 



P40.75 




0.65 



50 100 150 200 250 300 350 400 
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Fig. 2. The expected shape of R as function of Nw this time given directly by eq. 
© and data (SJ; no parameter p enters here. 



So far we have not yet used connection of R with EoS variables p and e as 
given by eq. Pj). According to it we can interpret the peculiar shape of curve 
obtained in Fig. ^ in terms of the type of EoS admitting it, in particular 
as due to a specific behaviour of p/e. Following therefore discussion in [2] 
we argue that this shape could be regarded as a signal of the existence of 
either the softest point of the corresponding EoS ^U] or its critical point 
jllj . This time it would show up as function of the number of participants 
(which can be translated into a volume of the reaction) rather than energy, 
as discussed in (2]. 

9 This could be regarded as a possible suggestions for experimentalists that everything 
needed for such discussion as presented here is measurable directly, without resorting 
to quantities like Ft or $, see also (17) for similar conclusions obtained at different 
circumstances. However, care must be excercised when following such approach be- 
cause usually bins in centrality are large and one can expect therefore contributions 
to fluctuations coming from fluctuations in read-offs of the ZDC calorimeter, which 
for a time being are not know. 
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Fig. 3. The same as in Fig. \I]p but this time using recent data on $ parameter 
obtained by CERES Collaboration J3| ■ The results are presented both as function 
of participants ((a)) and as function of energy ((6)). 



Actually, because of relation between Ft and measure $ of fluctuations 
|14j (see (|B.9|) in Appendix B), measurements of variable $ is equally good 
for the kind of analysis performed here. This is clearly seen in Fig. ^) and 
fusing NA49 data jU] and in Fig. where recent data on $ (cf., ^3]) were 
used as our input 10 . In Fig. |l]we show that data from clearly indicate 
that, as was already mentioned before, fluctuations of transverse momentum 
defined by variable u introduced in eq. (J7J) depend very weakly both on the 
energy and on the centrality of the collision (i.e., on the number of struck 
nucleons Nyy)- Therefore fluctuations of transverse momenta are practically 
irrelevant for the problem considered here, i.e., for the apparent structure 
seen in the EoS. Main effect is provided by fluctuations of multiplicity N, 
as presented by eq. ©• 

Let us close this Section with the following remark. Because fluctuations 
in J2pt are tightly connected with fluctuation of temperature T, therefore 
one can write that = and express parameter R in yet another form: 

Var(N) 

R = Var(T) ■ ( 15 ) 
T' 1 



10 



Notice that in both cases the lack of fluctuations in N, i.e., Var(N) = would 
immediately result in $ = Ft = 0. Lack of fluctuations means that also p — what 
(together with y — 0, cf. 11 H ) results in R — as well. 
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As shown in 20 11 relative fluctuations of temperature T are connected 
with the heath capacity and can be parameterized by the nonextensivity 
parameter q, 

Var{T) N w 



y2 



C v 



1, 



(16) 



which in our case is given by 



(N) 



0J + 



Var(N)' 



(17) 
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Fig. 4. Fluctuations of the transverse 
momentum as denned in eq. J7J as 
function of centrality and for different 
energies (data are from they are 
presented before removal of short range 
correlations). 



Fig. 5. Dependence of the NB [T3] 
multiplicity distribution parameter 1/k 
on the number of participants N\y for 
PHENIX 8 and CERES data. 



Notice that results shown in Fig. |1J i.e., that essentially uj = const, 
means that q depends only on the multiplicity (or on the number of struck 
nucleons), in fact q — 1 diminishes as 1/{N). Parameter q allows us then 
(by using eq. 0) to express parameter R in the form stressing its vital 
dependence on fluctuations of multiplicity, i.e., on Var(N)/(N): 

1 Var(N) 

R ~q-1 W ' ( } 



For other hints on nonextensivity in hadronic production processes and references to 
nonextensive statistics, see [21]. 
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As is well known these fluctuations lead to deviations from the poissonian 
form of multiplicity distributions in multiparticle production processes and 
result in broader distributions usually expressed by the so called Negative 
Binomial (NB) form and characterized by the parameter k 12 , which in our 
case is given by: 

1 Var(N) 1 . , 1 

Its dependence on N w for PHENIX data (at y/s = 200 GeV) and for 
CERES data (at \/s = 17 GeV) is shown in Fig. [5J As seen there param- 
eter 1/k decreases with increasing number of participants. It means that 
(cf. ^Hl) with increasing centrality fluctuations of the multiplicity become 
weaker and the respective multiplicity distributions approach poissonian 
form. Notice that because R < 1 the following condition, 

S + W <,_1, < 20 > 

must be satisfied. 



3. Nonmonotonicity and phase transition 

Let us now continue discussion of the parameter R from the point of 
view of its possible connection with the phase transition. Using known 
thermodynamical identities one can rewrite R (see Appendix C, all 
derivatives are for T = const and p = as) as 

( dlnV \ 2 
1 - Oyi EdWE 

Notice that because in the vicinity of critical point {dp/dV) T — > the 
parameter R reaches its maximum there: 

R^-mr. (22) 



(dEY 

[dVjr 



dE\ 2 \dlnEJ T 



Let us discuss eq. ()21(l in more detail concentrating on two cases important 
for us. First of all notice that for fixed energy y/s we have ^dE/E = dV/V = 
dN\y/N\y and with increasing iV^y the value of R decreases for £ = 1 and 

12 Cf. |19| for details. Connection between NB and nonextensivity in hadronic collisions 
represented by q > 1 has been discussed in [22]. 
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increases for £ = — 1. It means then that for fixed Nyy parameter R goes 
through its maximum when dlnV/dlnE changes sign. Actually exactly 
such change of sign is observed in compilation of heavy ion data on central 
collisions (for which Nyy remains approximately fixed) |15j . see Fig. HJ., 
from which one can deduce that 




d\nV J ci = -0.73 for s < s c , . 

m^E~{ c 2 = +0.15 for s>s c {Z6) 

Using now these values together with eq. (|2*Tj) we have obtained depen- 
dence of R on energy ^/s, cf. Fig. Eh- The parameters used were such 
that 

N w aT 2 2aT 



1 N W q-1 



30 



what corresponds to the reasonable values of a = 1/3, T = 0.14 GeV and 
q — 1 ~ 1/k = 0.003 13 . Notice that R shown there reaches its maximum 
for energy ~ 6 GeV. This " critical" value of energy is given by 

aT 2 c lC2 (M) 



q-1 N w ci + c 2 ' 



13 One should notice that in eq. (I2H we have total energy: E = ^fsNw /2. This is the 
origin of appearance of Nw in presenting R as function of ^/s. 
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where Cj denote the corresponding values of derivatives d In V/d In E as given 
in (|23j) . Actually, although the shape of -R in Fig. 03 resembles that ex- 
pected in [2], the height of its maximum is much smaller than expected 
there (to get the value observed in [2J one would have to shift it upward by 
~ 0.7). The only possible explanation we could offer at the moment is the 
observation that R displayed in Fig. 03 does not contain contribution from 
statistical fluctuations of multiplicity (which would give precisely the seek 
for value of 0.69 for the poissonian fluctuations). It would then mean that 
what we are calculating here are only changes in V . For V = const we are 
therefore getting R = 0. On top of that one has poissonian fluctuations, 
which for V =const provide, as mentioned above, the lacking ~ 0.7 on top 
of what is observed. In obtaining it we have increased fluctuations of mul- 
tiplicity N by poissonian component (adding to Var(N) the value of (JV) 
for constant volume V), i.e., we have used here R' instead of R as our input 
information, where 1/(1 — R') = 1/uj + 1/(1 — R). It is worth to notice that 
recent lattice calculations show behaviour of R similar to obtained by us 
here [TT] . 



4. Summary and conclusions 

We have shown that the recently proposed method of (almost) direct 
investigations of the EoS, especially finding a possible traces of phase tran- 
sitions to QGP phase of matter [2] by analysing energy dependence of some 
specific fluctuations observed in the multiparticle production data obtained 
in heavy ion collisions, can be further extended to include also fluctuations 
of different types than those proposed in |2j. Three examples were discussed 
here: (£) recent PHENIX Collaboration data [8] on F T ; (ii) recent NA49 
Collaboration data 0| on <3? and (iii) CERES Collaboration data |15l I13| 
on <£. In particular, we have derived formula (our eq. (jlOj) ) expressing pa- 
rameter R introduced in [2] by the measured quantity Ft (defined in Q 
and measured in [S]). Our results (see Fig. show that one observes 
similar characteristic structure in R as that expected in [2| but now being 
present at given fixed energy and for different centralities (expressed by the 
number of participants Nyy)- As is seen in Fig. Eland Fig. [I]p, the same 
type of analysis can be performed using as input the so called $ measure 
of fluctuations as provided by recent NA49 [S| and CERES data 13 . We 
would like to stress here that because, as is witnessed by results shown in 
Fig. 0J fluctuations in are rather irrelevant here, our results are almost 
entirely due to the fluctuations in multiplicity. It means therefore that they 
are not sensitive to flow phenomena. 

We have also established connection of R with fluctuations of temper- 
ature T described by nonextensivity parameter q [20]. The fact that data 
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|15j clearly indicate that in the measured centrality range fluctuations in 
transverse momentum are essentially constant, see Fig. 01 allows us to ex- 
press parameter R by combination of fluctuations of T (i.e., by parameter 
q) and fluctuations of multiplicity as given by the parameter 1/k of NB 
distribution [TJJI, see eq. ()19|) and Fig. |SJ 




10 2 5 10 2 2 5 10 3 2 5 10" 2 5 ' 10 2 5 10 2 2 5 lO* 2 5 lO* 2 5 

E = N w s 1/2 /2 [GeV] E = N w s 1/2 /2 [GeV] 



Fig. 7. EoS corresponding to results obtained in Fig. (solid line) and to those 
obtained from Fig. H] (points). Open points correspond to data from NA49 [HJ and 
full points to data from PHENIX [H]. See text for details. 

Finally, using some thermodynamical identities, we have expressed R in 
terms of parameter a = p/e defining form of EoS used and energy depen- 
dence of the reaction volume, cf. eq. (|21|). Using recent data on the later 
|15j (see Fig. [BJi) we have obtained the characteristic " shark fin" structure of 
the energy dependence of the parameter R shown in Fig. 03. It corresponds 
to the shape of EoS as given in Fig. UJ Such structure of R was predicted at 
[2] (albeit for ^/s replaced by the so called Fermi collision energy measure) 
as reaching value of R ~ 0.8 at its maximum. In our case the maximum is 
much lower. Actually, what is shown in Fig. [JJ are two situations, which 
can emerge from definition of R as given in eq. ifTJl. Namely, as one can 
notice variation of R can originate either in variations of p/e or in variations 
dV/dT (in extremal cases). For dV/dT = we are then getting from data 
p/e, and this is the case of Fig. [7^. On the other hand, for p/e = 1/3, we 
are getting dV/dT, as in the Fig. [JJa. 

We would like to close with the following remark. Our result presented 
in Fig. [JJ shows not one but two "soft points" or "critical points" of EoS 
(so far we cannot distinguish here between these two possibilities), located 
at different energies: one at y/s ~ 6 GeV (obtained from data on V at 
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different energies |15| and from data on $ for different centralities [0]) and 
another one at yfs ~ 100 GeV (obtained from data on fluctuations mea- 
sured by quantity Ft for different centralities |5j). We cannot so far offer 
explanation of this result. Interpreting it, however, in the original spirit 
of looking for the possible signals of QGP (but keeping in mind our reser- 
vation concerning this point expressed above) we could say that it looks 
like, with increasing energy, first occurs a kind of QGP composed mostly 
of dressed quarks, which at higher energies is followed by a true QGP con- 
taining also liberated gluons 14 . What we see here most probably indicates 
some additional change in energy dependence of V presented in Fig. HJ., 
which would take place at \fs ~ 100 GeV. If true it would indicate that V 
is not increasing with energy since that point, on the contrary, it probably 
decreases a little (what seems to be confirmed by the first data from RHIC 
at 200 GeV [2SD- One should also keep in mind, before further speculations, 
that the assumed here access to the EoS by use of fluctuations in hadronic 
production neglects altogether the possible evolution processes between the 
QGP and state of freely streaming hadrons observed experimentally. This 
remark applies also to all works of this kind, like [2]. It means therefore 
that final conclusions could only be drawn when other approaches, as for 
example that discussed in [7j, are also confronted with all available data, as 
has been done here. 



Appendix A 

We shall, for completeness, derive here formula (Q). We shall consider 
in what follows our hadronizing system as statistical system in equilibrium 
without specifying its details. Thermo dynamical potential of such system 
is 

ft = E + pV - TS = fiN. (A.l) 

It is then connected with chemical potential n with N denoting number of 
particles in the system under consideration. For system in equilibrium but 
with varying number of particles N is given by condition dVt/dN = or, 
equivalently, [i = 0. For conditions of statistical equilibrium we have then 
(here e and s are the corresponding densities of energy and entropy) 

£ + p - Ts = 0. (A.2) 

Using (|A.2|) together with thermodynamic identity 

dE = -pdV + TS (A.3) 

14 Notice that data from |]T5| indicate that volume of interaction grow much slower than 
linearly with energy, i.e., the energy densities obtained at different energies increase 
substantially: e ~ 0.1s ' 41 . 
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we get 

de + dp = Tds + sdT. (A.4) 
Denoting a = Cq = dp/de = p/e we obtain from (jA.2|) and l|A.4|) 

ds 1 dT , s . . 

— = and de = -dT. (A.5) 

s a T a 

Because from (jA.2j) we get e = Ts — p = Ta(de / dT) — ae, then 

* = I±^. (A.6) 
e a T 

From equations (|A.5|) and (|A.6() one gets that ratio of fluctuations of den- 
sities is given by 

jds/s) 2 1_ 

de/e) 2 ~ (1 + a) : 

Because 5 = and E = eV and, respectively, dS = dsV + sdV and 
dE = Vde + edV, then 



d5 1 dT / IW\ d£ 1 + aeff/ a T dV 

- = a^{ 1 + a vdr) and ^ = ^r^( v 1 + TT^Fdr;' 

(A.8) 

and the corresponding ratio of fluctuations of entropy and energy is equal 

(dS/S) 2 ( a \~ 2 TdV dlnV 

R = i- L '— =1-| where 



(dE/E) 2 \ l + a</)J V dT dlnT 

(A.9) 

For dV = one gets result obtained already in 2 j. However, in general 
</> 7^ and fluctuations of the interaction volume are also important. Eq. 
jl can be also rewritten as 



- = -4> + -^=, (A.10) 
a 1 - VR 



which connects EoS (represented by a) with changes in the reaction volume 
and with fluctuations described by the parameter R. 

Appendix B 

We shall derive here eq.fjHJ). Let us first notice that Ft defined in @ is 
related to other measure of fluctuations, the so called ^-measure introduced 
in [Hj and used in piT3]. 



(Z 2 ) 
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For a measured quantity x (here identified with transverse momentum of 
produced particles, Xj = pxi) one has: 

z = x — (x), i^ 2 ^) = Var(x) = (x 2 ^j — (x) 2 
N N 

z = 52zi = J2 x i- N ( x >> ( B - 2 ) 

i=l i=l 

( Z2 ) = ((l>) )-2(x)^£^ + (iV 2 )(x) 2 . 

To get this result one uses following reasoning. If variables x and N are 
described by distributions characterized by the respective generating func- 
tions f{t) and h{t) then and variable N by generating function h(t) then 
variable £ = Yli=i x i ls described by generating function G(t) = h[f(t)]. It 
means therefore that 

(Z)=G'(l) = (N)(x) (B.3) 

and 

(£ 2 ) = G"(l) + G'(l) = (N 2 )(x) 2 + (N)Var(x) (B.4) 

leading to 

Var{£) = (N)Var(x) + (x) 2 Var(N). (B.5) 

To characterize correlations between variables N and £ one has to introduce 
a correlation coefficient p G [—1, 1], which will be, in what follows, our free 
parameter (our minimal input). One can therefore write: 



N \ IN 



jvY>= (AO (!>+/> 



\ 



N 



Var{N)Var[J2xi). (B.6) 



\ i=l / \i=l / 

The last term in (IB. 61) can be written as 



\i=l 



\ 



N 



Var(N)Var [J2 x *j = V ' (N)Var(N)Var(x) + (x) 2 Var 2 (N). (B.7) 
Substituting now (|B.6|) to (|B.2|) and making use of eq. (|B.3I) - (|B,5|) one gets 
{Z 2 ) = {N)Var(x) + +2{x) 2 Var(N) 



I (N) Var(x) 
P VVar(N) (x) 2 



(B. 



Substituting this to eq. (|B.1|) and making use of the relation between Ft 
and <!> derived in ^f)j, namely that 
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one gets 
F T = -1 + 



1 + 2- 



(x) 2 Var(N) 



\ * ' "Far(a;) (N) 
which is eq.® we were looking for. 



/ (N) Var(x) 
P V Var(N) (x) 2 + 



(B.10) 



Appendix C 

Let us start with change in the energy E of the system, which can be 
written as (Cy is corresponding heat capacity) 

dp ' 



AE 



(™) AV + (™) AT 
\dVJ T \9TJv 



AV + CyAT. (C.l) 



Squaring it and averaging while remembering that [2H fluctuations of the 
volume and temperature are given by, respectively, 

'dV\ T 2 

I a.nrl ( ATI 2 = 

C v 



{AVf 



- Tl ^) T and ^ 



(C.2) 



and are statistically independent, i.e., AT6V = 0, one gets 



Var(E) = (AE 2 ) 



P 



(C.3) 



Similarly, fluctuations of number of particles can be described by formula 
From (|C.3|) and (|C.4j) one gets 



(C.4) 



R 



Var(N) 
Var(E) 



CyT 

~7T- 



Because 



then (from now on all derivatives are for T = const) 



(C.5) 



(C.6) 
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Because of EoS one has that 



and 



dp dp dV a 



(C.8) 



p = ae 



dE dV dE V 



what immediately leads to eq. (|21|k 



R = 




(C.9) 



i /~~t rp a dV 

1 - c v i vm 
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